Chiral transition in a strong magnetic background 
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The presence of a strong magnetic background can modify the nature and the dynamics of the 
chiral phase transition at finite temperature. We compute the modified effective potential in the 
linear sigma model with quarks to one loop in the MS scheme for Nf = 2. For fields eB ^ 5m^ 
and larger a crossover is turned into a weak first-order transition. We discuss possible implications 
for non-central heavy ion collisions at RHIC and LIfC, and for the primordial QCD transition. 
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I. INTRODUCTION 

Strong magnetic fields can produce remarkable phys- 
ical effects. One of the most spectacular examples is 
provided by magnetars |JL|. More recently, it has been 
proposed that strong magnetic fields could play a very 
important role in the physics of high-energy heavy ion 
collisions, affecting observables in the case of non-central 
collisions at the Relativistic Heavy fon Collider (RHfC), 
at Brookhaven, and the Large Hadron Collider (LHC), 
at CERN, and providing a possible signature of the pres- 
ence of CP-odd domains in the presumably formed quark- 
gluon plasma (QGP) phase [2]- In that case, one would 
reach magnetic fields B ~ 10^^ Gauss, which correspond 
to eB ~ 6to^, where e is the fundamental charge and 
is the pion mass. This gives a very intense magnetic 
background from the point of view of quantum chromo- 
dynamics (QCD) scales (iTj . 

In this paper we investigate the effects of a strong 
magnetic background on the nature and dynamics of 
the chiral phase transition at finite temperature, T, 
and vanishing chemical potential. As a framework, we 
adopt the hnear sigma model coupled to quarks (LSM^) 
with two flavors, Nj = 2 This effective theory 

has been widely used to describe different aspects of 
the chiral transition, such as thermodynamic properties 
i, i, i, 0, i, i, 113, im, E M,M and the nonequilib- 
rium phase conversion process [l5|, as well as combined 
to other models in order to include effects from confine- 
ment 

In the limit of strong magnetic fields, we find that the 
nature of the chiral transition is modified. The origi- 
nal LSMg, with no magnetic field, yields a crossover at 
T = Tc ^ 150 MeV for a choice of coupling to quarks that 
provides a reasonable mass for the constituent quarks at 
zero temperature 8]. The presence of a strong mag- 
netic background turns this picture into a weak first- 
order phase transition. As a consequence, the dynamics 
that follows a rapid supercooling, such as the one that 
presumably happens in a high-energy heavy ion collision, 
can be dramatically modified due to the presence of a 
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barrier in the effective potential. Even if the barrier is 
quite small, as in the case of a weak first-order transi- 
tion, its effect can be very significant for the dynamics, 
holding the system in the false vacuum until it reaches 
the spinodal instability and explodes [Tl ]. 

From the theoretical point of view, the non-trivial role 
played by magnetic fields in the nature of phase transi- 
tions has been known for a long time [l^ . Modifications 
in the vacuum of quantum electrodynamics (QED) and 
QCD have also been investigated within different frame- 
works, mainly using effective models [20|, l2l|, [H, [H, [13, 
m, , especially the NJL model [13] , and chiral pertur- 
bation theory [H, [1^, [s^ , but also resorting to the quark 
model [sill and certain limits of QCD [32*1 . Most treat- 
ments have been concerned with vacuum modifications 
by the magnetic field, though medium effects were con- 
sidered in a few cases. More recently, interesting phases 
in dense systems [s^ , as well as effects on the dynamical 
quark mass [3^ and on the quark-hadron transition (35| 
were also considered. 

The paper is organized as follows. Section II presents 
briefly the low-energy effective model adopted in this pa- 
per, as well as a discussion of the approximations re- 
quired. In Section III we show our results for the modified 
effective potential, in the presence of a strong magnetic 
background. Phenomcnological consequences are illus- 
trated and discussed in Section IV. Section V contains 
our conclusions and outlook. 



II. EFFECTIVE THEORY 

To describe the chiral phase structure of strong inter- 
actions at finite temperature we adopt the LSM^ , defined 
by the following lagrangian 

C ^ tpf [ij'^d^ - .g(o- + il5T ■ 7?)] ipf 

+ l{d^ad^a + df,TTd''n)-V{a,TT) , (1) 

where 

V{a,Tf) = ^{a^+Tr^-vY-ha (2) 

is the self-interaction potential for the mesons, exhibiting 
both spontaneous and explicit breaking of chiral symme- 
try. The Nf = 2 massive fermion fields ipf represent the 
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up and down constituent-quark fields ip = {u,d). The 
scalar field a plays the role of an approximate order pa- 
rameter for the chiral transition, being an exact order 
parameter for massless quarks and pions. The latter are 
represented by the pseudoscalar field tt ~ (tt^, 7r+, 7r~), 
and it is common to group together these meson fields 
into an 0(4) chiral field = (cr, tt). 

In what follows, we implement a simple mean- field 
treatment with the customary simplifying assumptions 
(see, e.g., Ref. @). Quarks constitute a thermalized 
fluid that provides a background in which the long wave- 
length modes of the chiral condensate evolve. Hence, at 
T = 0, the model reproduces results from the usual LSM 
without quarks or from chiral perturbation theory for the 
broken phase vacuum [36| . In this phase, quark degrees of 
freedom are absent (excited only for T > 0). The a field 
is heavy, ~ 600 MeV, and treated classically. On the 
other hand, pions are light, and fluctuations in tt"*" and 
TT~ couple to the magnetic field, B, as will be discussed 
in the next section, whereas fluctuations in tt^ give a B- 
independent contribution that we ignore, for simplicity. 
For T > 0, quarks are relevant (fast) degrees of free- 
dom and chiral symmetry is approximately restored in 
the plasma for high enough T. In this case, we incorpo- 
rate quark thermal fluctuations in the effective potential 
for cr, i.e. we integrate over quarks to one loop. Pions 
become rapidly heavy only after Tc and their fluctuations 
can, in principle, matter since they couple to B. How- 
ever, they will play a minor role in this case, as will be 
clear later, so that we ignore their thermal fluctuations 
in this section, for simplicity. 

The parameters of the lagrangian are chosen such that 
the effective model reproduces correctly the phenomenol- 
ogy of QCD at low energies and in the vacuum, in the 
absence of a magnetic field, such as the spontaneous 
(and small explicit) breaking of chiral symmetry and ex- 
perimentally measured meson masses. So, one has to 
impose that the chiral SUl{2) (g) SUr{2) symmetry is 
spontaneously broken in the vacuum, and the expecta- 
tion values of the condensates are given by (cr) = /^r 
and (tt) = 0, where /^r = 93 MeV is the pion decay 
constant. The explicit symmetry breaking term is de- 
termined by the PC AC relation which gives h — /tt^Jt-^, 
where m^r ~ 138 MeV is the pion mass. This yields 
v"^ = fl ~ m^/X. The value of A = 20 leads to a a- 
mass, = 2A/^ + "^^1 equal to 600 MeV. In mean field 
theory, the purely bosonic part of this Lagrangian ex- 
hibits a second-order phase transition [33| at Tc = \f2v 
if the explicit symmetry breaking term, /ig, is dropped. 
For hq 7^ 0, the transition becomes a smooth crossover 
from the restored to broken symmetry phases. For g > 0, 
one has to include a finite-temperature one-loop contri- 
bution from the quark fermionic determinant to the ef- 
fective potential. The choice g — 3.3, which yields a rea- 
sonable mass for the constituent quarks in the vacuum, 
Mq = g/jr, leads to a crossover at finite temperature and 
vanishing chemical potential, and will be adopted in this 
paper |42|]. 



Standard integration over the fermionic degrees of free- 
dom [36j . using a classical approximation for the chiral 
field, gives the following formal expression for the effec- 
tive potential in the a direction: 



K//(T,a) = F(a) 



T 
V 



In det 



(G^i+M,(a)) 
T 



(3) 



where Vip) is the classical self-interaction potential for 
the mesonic sector in the cr direction, Ge is the fermionic 
Euclidean propagator, Mq is the effective fermion mass 
in the presence of the chiral field background, and V is 
the volume of the system. To one loop, one has Veff — 
y((/)) -t- Vq(0), where the contribution from the quarks 
reads 



~VqT 



(2^ 



In 



(4) 



where Uq = 24 is the color-spin-isospin-baryon charge de- 
generacy factor, Ek{a) = {P + M^{a)y/^, and M{a) = 
g\(j\ plays the role of an effective mass for the quarks. 
The net effect of this term is correcting the potential for 
the chiral field, approximately restoring chiral symmetry 
for a critical temperature Tc ^ 150 MeV Q. The effec- 
tive potential Vef / for several values of the temperature 
is displayed in Fig. 1, assuming g — 3.3. 
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FIG. 1: Vsff{a) for different values of the temperature for 
g = 3.3, in the absence of a magnetic field. 



III. MODIFIED EFFECTIVE POTENTIAL 

Let us now assume that the system is in the presence of 
a strong magnetic field background that is constant and 
homogeneous. For definiteness, let us take the direction 
of the magnetic field as the ^-direction, B = Bz. One can 
compute the modified effective potential by redefining 
the dispersion relations of the scalar and spinor fields in 
the presence of B, using the minimal coupling shift in 
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the gradient and the field equations of motion. For this 
purpose, it is convenient to choose the gauge such that 
A^^^{A\A)^{0,-By,Q,Q). 

For scalar fields with electric charge g, one has 



{d^ +m^)4> = 0, 



(5) 
(6) 



After decomposing cf) into Fourier modes, except for the 
dependence in the coordinate y, one obtains 



2m 

2m 



? 9 ? 

2m 



Px 

qB 



^(y) = o. 



(7) 



which has the form of a Schrodinger equation for a har- 
monic oscillator. Its eigenmodes correspond to the well- 
known Landau levels 



2 2 
POn ~ Pz 



2m 



(8) 



where lob = \q\B /m and n is an integer, and provide the 
new dispersion relation: 



Pon ^ pi 



{2n+\)\q\B. 



(9) 



One can proceed in an analogous way for fermions with 
charge q. From the free Dirac equation {i^^d^ — m)il} = 
0, and the shift in 9^, one arrives at the following 
Schrodinger equation 



2m 

q^B^ 
2m 



9 2 



m 



qBa 



2m 

2' 



y 



Px \ 
qB) 



u,(y) = 0, (10) 



which yields the new dispersion relation for quarks: 

1.^ + {2n ^ \ - a)\q\B . (11) 



-POn 



Pz 



It is also straightforward to show that integrals over 
four momenta and thermal sum-integrals acquire the fol- 
lowing forms, respectively: 





^ \q\B ^ r dfco 
27r ^ / 27r 

n— 


dk^ 




27r ' 


d^k 

(27r)3 




f dk, 
27r 



(12) 
(13) 



where n represents the different Landau levels and £ 
stands for the Matsubara frequency indices [s^ . 

According to the assumptions of our effective model, 
the vacuum potential will be modified by a contribution 



from the charged pions that couple to the magnetic field. 
To one loop, this is given by the following integral: 



n=0 '' 



1/2 



(14) 



For sufficiently high magnetic fields, eB >> mj, and 
ignoring contributions independent of the chiral conden- 
sate, this simplifies to the following renormalized result 
0: 

e+e=-^log2, (15) 
where we have added the contributions from 7r+ and 7r~. 
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FIG. 2: Veff{a,T,B) for T 
magnetic field. 



and different values of the 



Results for the modified vacuum potential are dis- 
played in Fig. 2 for different values of the background 
magnetic field. As one can see from the plot, the pres- 
ence of the magnetic field enhances the value of the chiral 
condensate and the depth of the broken phase minimum 
of the modified effective potential, a result that is in line 
with those found within different approaches (see, for in- 
stance, Refs. [m, 1131 )• This could be an indication 
that the presence of a strong magnetic background might 
favor a weak first-order transition over the crossover that 
is obtained for _B = 0. To investigate this point, one has 
to compute the modified thermal contributions to the 
effective potential and study the process of chiral sym- 
metry restoration. 

Charged pions contribute to the one-loop thermal cor- 
rections via 



\q\BTf. 

271 ^ 



dk 
2^ 



log(l 



(16) 



whereas the contribution from quarks reads 



71=0 



dk, 

2^ 



21og(l + e 



(17) 



where, for quarks, one still has to sum over colors, which 
brings a factor Nc, and over flavors. The latter sum, 
differently from what happened in the previous section, 
does not bring and overall Nf factor, since the electric 
charges are not equal, and their absolute values add up 
to one. In the presence of a nonzero B, spins are also 
treated separately, as is clear in the expression above. 
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FIG. 3: Veff{(J,T,B) for eB — 5mJ and different values of 
the temperature. 

For sufficiently high magnetic fields, not only compared 
to the pion mass but also to the effective quark mass and 
to the temperature, and again ignoring contributions in- 
dependent of the chiral condensate, we obtain for charged 
pions and quarks, respectively: 



'(27r)3/2 



(&) 


n 


ml ' 









(18) 



eBT' 



+ 00 



dx log I 



ieB 



(19) 



One can notice that the pion thermal contribution 
and part of the quark thermal contribution are exponen- 
tially suppressed for high magnetic fields, which essen- 
tially comes from an increase in the effective mass due to 
the magnetic field - e-^^-f^/"^ ~ ^-#V^/t^ Therefore, 
within these approximations, the total modified effective 



FIG. 4: Veffia,T,B) for eB 
the temperature. 



potential is given by 



10m„ and different values of 



27V, 



eBT^ 



dx log I 



_ log2 

e-V^""^^] . (20) 



The modified effective potential, Veff{a, T, B), is plot- 
ted for different temperatures for eB — 5m^, eB — lOmJ 
and eB = 20m^ in Fig. 3, Fig. 4 and Fig. 5, respec- 
tively. For eB = 5rn^, one verifies that the critical tem- 
perature becomes higher, > 200 MeV. Moreover, a 
closer look uncovers the presence of a tiny barrier, sig- 
naling a first-order phase transition. A zoom for a couple 
of temperatures close to Tc is shown in Fig. 6. In the 
case of higher magnetic fields, the existence of the bar- 
rier is rather clear, showing that the presence of a strong 
magnetic background can really modify the nature of the 
chiral transition, in this case from a crossover to a first- 
order transition whose strength depends on how intense 
B is. In Figs. 4, 5 and 6, we considered the same set of 
temperatures to allow for a more direct comparison. 



Another important feature that comes out of this com- 
parison is the fact that for eB — lOm^ the critical tem- 
perature drops again (Tc < 140 MeV) as compared to the 
case in which eB — 5m J. This phenomenon continues to 
happen for higher B (see Fig. 6), showing that there is 
a non-trivial balance between temperature and magnetic 
field that can bring some richness to the chiral T — B 
phase diagram that should be explored. The presence 
of the barrier is not qualitatively modified, though, and 
the "conversion" to a first-order transition seems to be 
robust. 
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FIG. 5: Veff(a,T,B) for eB = 20ml and different values of 
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FIG. 6: Veffip-, T, B) for eB = 5m^ and different values of the 
temperature (zoom near the critical temperature to exhibit 
the barrier). 



IV. PHENOMENOLOGICAL CONSEQUENCES 

Besides the pure theoretical interest in the modifica- 
tion of the nature of the chiral phase transition by the 
presence of a strong magnetic background, one should 
consider possible phenomenological implications within 
the reach of present experiments. At RHIC, estimates by 
Kharzeev, McLerran and Warringa give eB ^ 5.3rri^. 
For the LHC, there is a factor Zpi,/Zau — 82/79 coming 
from the difference in charge, and some small increase in 
the maximum value of eB due to the higher center-of- 
mass energy (as observed for RHIC So, it is reason- 
able to consider eB 6m^ as a representative estimate. 

We show in Fig. 7 a zoom of the effective potential for 
eB 6mJ for a temperature slightly below the critical 
one, whereas the general picture is not appreciably differ- 
ent from the case eB ^ 5mJ, and can be seen in Fig. 3. 
Comparing Figs. 1 and 3, one can establish the basic dif- 
ferences brought about a magnetic field of the magnitude 
that could possibly be found in non-central high-energy 
heavy ion collisions. One moves from a crossover sce- 



nario to that of a weak first-order chiral transition, with 
a critical temperature ~ 30% higher. 

In the case of a crossover, one expects the system to be 
smoothly drained to the true vacuum, with no formation 
of bubbles or spinodal decomposition. In the case of the 
first order transition depicted in Fig. 3, despite the fact 
that the barrier is quite small, part of the system will be 
kept in the false vacuum for a while, a few bubbles will be 
formed, depending on the intensity of the supercooling, 
and spinodal instability will then set in [Tl ]. 

The case of the LSM, with a coupling constant g = 5.5 
exhibits a first-order phase transition even for B — 0, 
and has been often used to study nucleation and spin- 
odal decomposition phenomena in the chiral transition 
S i, 0, [11 Q E [i3- There, one finds a barrier 
~ 0.25 in the effective potential for temperatures close 
to Tc, in units of T^, and it has been shown that the 
system stays mostly apprisionated in the false vacuum 
until it reaches the spinodal explosive phase conversion 
regime joj. For g — 3.3, we saw that there is a very weak 
first-order chiral transition for i? > 0, with a barrier an 
order of magnitude smaller (see Fig. 7). Nevertheless, 
even such small barriers can hold the system in the false 
vacuum until the spinodal instability in the case of a fast 
supercooling, as expected for high-energy heavy ion col- 
lisions Ts*]. Therefore, the existence of a strong magnetic 
background can bring significant effects to the dynamics 
of the chiral transition. 
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FIG. 7: Veffia, T, B) for eB = 6ml for T = 180 MeV. 



V. CONCLUSIONS AND OUTLOOK 

Lattice QCD indicates a crossover instead of a first- 
order chiral transition at finite temperature and vanish- 
ing chemical potential. However, a strong magnetic back- 
ground might modify this situation. Our computation of 
the modified effective potential clearly shows the forma- 
tion of a barrier that seems to be robust, in the case of 
magnetic fields that are larger than the other relevant 
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scales, namely the pion mass and the temperature. On 
the other hand, it is clear that further investigation of the 
low magnetic field regime at finite temperature, for both 
B < T and i3 ~ T, in the context of the chiral transition 
is called for, since the phase diagram seems to exhibit 
non-trivial features in this region of parameters which is 
interesting for critical phenomena in strong interactions. 

For heavy ion collisions at RHIC and LHC, the barrier 
in the effective potential seems to be quite small. Nev- 
ertheless, it can probably hold most of the system in a 
metastable state down to the spinodal explosion, imply- 
ing a different dynamics of phase conversion as compared 
to the crossover scenario. Moreover, the phenomenol- 
ogy resulting from varying T and B seems to be rich: 
there seems to be a competition between strengthening 
the chiral symmetry breaking via vacuum effects and its 
restoration by the thermal (magnetic) bath. In particu- 
lar, non-central heavy ion collisions might show features 
of a first-order transition when contrasted to central col- 
lisions. 

The dynamics of the chiral phase transition will cer- 
tainly be affected by the presence of a strong magnetic 
background, since the nature of the transition is modi- 
fied. In order to study its impact on the relevant time 
scales of this process, one has to perform real-time sim- 
ulations of the evolution of the order parameter in the 
presence of the modified effective potential, e.g., in a 
Langcvin framework. Results in this direction will be 
presented elsewhere [38| . 

Of course, the description presented in this paper is ad- 
mittedly very simple. In actual heavy ion collisions, the 



magnetic field is neither uniform nor constant in time, 
and the values quoted here represent estimates of its max- 
imum magnitude. In fact, the field intensity rapidly de- 
creases with time, and any of the effects discussed here 
should manifest in the early-time dynamics. Moreover, 
if there is a a fast variation of B with time, there will 
be an induced electric field that could break the vacuum 
through the Schwinger mechanism of pai r production, 
and this would modify the condensate [30| and the vac- 
uum effective potential. Therefore a more realistic de- 
scription of the external electromagnetic field effects on 
the chiral transition must be pursued before one is able 
to make explicit phenomenological predictions. On the 
other hand, we believe that the results presented here, 
even if seen as just qualitative, are encouraging. In par- 
ticular, the possibility that strong magnetic fields can 
turn a crossover into a first-order transition might be 
relevant for the physics of the primordial QCD transi- 
tion, eventually reconciling lattice results that predict a 
crossover [s^ to the description of phase conversion via 
bubble nucleation that assume a first-order phase transi- 
tion [13 . 
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